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Abstract 

The two B decay processes B° — ► D*~ir + and B° — > D*~£ + vi have been studied by means of a 
partial reconstruction technique using a data sample collected with the BABAR detector at the 
PEP-II storage ring. To increase statistics, only the soft ir~ from the decay D*~ — > ir~D° was 
used in association with either an oppositely-charged high-momentum pion or lepton. Events were 
then identified by exploiting the constraints from the simple kinematics of T(4S) decays. A clear 
signature is obtained in each case. The position of the B° decay point was obtained from the 
reconstructed -K + (£ + )7r~ vertex. The position of the other B° in the event was also determined. 
Taking advantage of the boost given to the T(4S) system by the asymmetric beam energies of 
PEP-II, the lifetime of the B° meson has been measured from the separation distance between the 
two vertices along the beam direction. The preliminary results are: 

t b o = 1.55 ± 0.05 ± 0.07 ps, 
t b o = 1.62 ±0.02 ±0.09 ps, 

respectively for the B° — > D*~ir + and B° — > D*~£ + V£ channels. 
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1 Introduction 



This document presents two analyses designed to select large samples of B mesons with inclu- 
sive reconstruction techniques. The first method finds B° — > D*~it + events and the second 
B° — > D*~£ + vi; while the two techniques are different in detail, they both share the common 
feature of making no attempt to reconstruct the D° produced in the D*~ — > D°tt~ decay, thereby 
achieving high efficiency comparing to the exclusive reconstruction. 

Taking advantage of the boost of the T(AS) system at PEP-II, the lifetime of the B° meson 
is determined separately, with a good statistical precision, using both samples of events. Tagging 
the flavor of the second B° in the event will eventually also allow a measurement of the mixing 
parameter Am^. 

1.1 B° — > D*~ir + inclusive reconstruction 

The decay B° — > D*~ir + is interesting for many reasons, but particularly since it may exhibit CP 
violation [|]|. Since the asymmetry is expected to be small, this is a long term goal and the first 
steps toward it are the measurement of the B® lifetime and mixing using this process. 

The full reconstruction of the decay chain B° — ► D*~n + , D*~ — » D°tt~, where the D decays 
in the mode D° — » K~tt + , K~tt + tt°, K~tt~tt + tt + , K°ir~TT + , allows one to evaluate detector 
performance and to extract the signal with an excellent signal over background ratio. However the 
loss in statistics is substantial and may limit a CP violation measurement given the small expected 
asymmetry. 

An alternative approach exists for extracting the signal without reconstructing the decay of 
the D° meson. This partial reconstruction method has already been used successfully by other 
experiments [Q, || and allows an increase of the size of the reconstructed sample by about an order 
of magnitude while maintaining the background at a reasonable level. 

No attempt is made to reconstruct the D° decays. Therefore, one searches for a pair of 
oppositely-charged pions (jTf, tt s ) and, assuming that their origin is a B° meson, calculates the 
missing invariant mass which should be the D° mass if our hypothesis was correct. Without the 
constraint of the D° mass, the direction of the B meson is unknown. Although its angle with 
respect to 7rj direction can be deduced, the angle (ft around this direction is undetermined. For 
more details see [§, ||. 

Using the beam energy constraint, the missing mass, M m i ss , is computed from the energy and 
momenta of the two reconstructed pions. Since this still depends on the unknown angle (ft of the 
B° momentum, in this analysis, M m i ss is defined as the average of the maximum and minimum 
value of M m i ss over all (ft angles. 

A Monte Carlo study shows that the resolution on M m i ss obtained in this way is of the order of 
3MeV/c 2 , and that it is dominated by the tracking precision for the measurement of the two pions. 

1.2 B° — > D*~£ + ui inclusive reconstruction 

The B° semileptonic decay process B° — > D*~£ + U£ is characterized by high branching ratio and 
clear experimental signature. These facts allow the selection of large event samples with small 
background contamination. Therefore this process has been widely used in the past to determine 
several B° properties, along which its lifetime, r^o , and the B°B° mass splitting, Am^ [0]. In 
those measurements the D*~ particle was reconstructed through its decay D*~ — > D°ir~, and 
the D° was identified by means of exclusive reconstruction of a few final states which provide high 
selection efficiency and high background rejection. 
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Due to the restricted phase space of the D*~ — > D°ir~ transition it is however possible to tag 
this process by reconstructing only the low energy charged pion (hereafter called ir* ) from the 
D*~ decay, without reconstructing the D°. This "inclusive" approach allows us to select an 
event sample about ten times more abundant than the exclusive analysis. This method has been 
extensively used in the past by the ARGUS and CLEO collaborations at the T(45) to determine 
several B° and D° properties §, §, and by DELPHI and OPAL at LEP to measure \V c b\ and the 
B° — > D*~e + vi branching fraction |?], @]. Measurements of the B° lifetime and oscillation frequency 
with this approach have only been performed by the DELPHI collaboration at LEP [||, [l^] and are, 
at present time, among the most precise determinations of those quantities. 

2 Detector and data 

The data used in this analysis were collected with the BABAR detector (for a detailed description see 
[|TT[| ) at the PEP-II storage ring. The statistics analyzed corresponds to an integrated luminosity of 
7.9 (7.1) fb _1 of data collected on the T(45)peak, and 1.2 (1.0) fb _1 collected below peak, for the 
B° — * D*~ir + (J3 — ► D*~l + vp) analysis. The off-resonance data were used for background study 
purposes. 

Samples of about one million B°B° and B + B~ simulated events were also analyzed; in addition, 
about 28000 (175000) B°B° were generated, were one of the B° mesons was forced to decay to 
B° -> D*-tt + (D*-£+u). 

3 B° to D*~7r + analysis and results 

3.1 Inclusive reconstruction of B° — ► D*~ir + 
3.1.1 Track and signal selection 

The fast and slow pion have been selected by requiring Drift Chamber (DCH) and Silicon Vertex 
Detector (SVT) quality cuts, respectively. In addition, the fast pion should have a momentum in 
the center of mass frame between 2.114 and 2.404 GeV/c, and should not to be identified as a muon 
or as an electron. The slow pion should have a momentum in the laboratory of at least 50 MeV/c. 
A further set of cuts were applied to reject electrons and kaons from the fast and slow pion tracks 
using the calorimeter and the Detection of Internally Reflected Cherenkov light (DIRC). 

In order to obtain the signal, pairs of oppositely-charged tracks were required and the inclusive 
reconstruction method was applied, leading to a value for M m i ss . An optimized selection criterium 
is necessary in order to reduce the combinatorial background, mainly from continuum events. 

In order of decreasing discriminating power, the quantities used for signal selection are: 

1. i?2) the normalized second Fox- Wolfram moment, is required to be less than 0.35. 

2. Requiring that no other tracks be in a cone of opening angle 0.4 rad centered on the fast pion 
candidate momentum vector in the T(4S) system. 

3. Two sets of particles are defined: all tracks and calorimeter clusters, excluding the fast and 
slow pion (SET1), and a set where all particles in a cone around the reference D° direction 
are also excluded (SET2) (see section |3.2.2| ). A single discriminating variable is constructed 
using the Fisher method, based on 15 input variables. (For SET1, the scalar sum of the 
CM momenta of tracks and showers in nine 20° angular bins around the fast pion, and 
the sphericity. For SET2, the angle of the sphericity axis and the angle of the particle of 



8 



highest energy with respect to the fast pion, and the combined mass and momentum of all 
the particles.) 

4. The cosine of the helicity angleQ is required to be larger than 0.4 in absolute value. 

The efficiency for these cuts as determined using Monte Carlo simulation is 27.5%. 

The distribution of the missing mass obtained on the data, applying only the R2 cut, is shown 
in Fig. [j]. The fitted signal is 10992 ± 235 events. Applying all the cuts (see also Section 6.4) the 
distribution in Fig. is obtained, containing 1696 ± 107 events. 
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Figure 1: Missing mass distribution for par- 
tially reconstructed B° — > D*~ir + events 
from data. Here only the cut on R2 is ap- 
plied. 



Figure 2: Missing mass distribution for par- 
tially reconstructed B° — > D*~ir + events 
from data after all the cuts are applied. 



No peak is observed in the corresponding distribution for wrong charge combinations, nor in 
the off-resonance data. 



3.1.2 Signal and background composition 

Other B decay modes may produce an enhancement at the end of the missing mass spectrum in 
the same region where the signal is expected. An obvious example is the decay B° — > D*~ p + . 

To study this effect, we have fit separately the distribution obtained for generic B mesons and 
charged B mesons from Monte Carlo simulation. Contributions to the missing mass spectrum come 
from a) B° — > D*~ir + signal, b) combinations of a true soft pion from a D with a pion coming 
from a B meson or from a direct p, c) combinations of a true soft pion with a direct lepton, and 
d) all other combinations. Only the categories a), b) and c) show clearly a peak at the end of 
the spectrum. The fit gives 461 ± 26 events, which is consistent with the sum of the contribution 
from a) (365 events) and b) (100 events). Therefore, we estimate that the signal contribution is 
enhanced by 26% from other B° decay modes. 

No clear peaking is observed for charged B mesons. The estimate for their fraction in the signal 
sample is 6.0 ± 3.1%. 

1 angle between the soft pion and the D* flight direction in the rest frame of the latter 
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Figure |3| shows a comparison between data and Monte Carlo distributions for missing mass, 
after all cuts. The Monte Carlo is normalized to the data. 

Table 1: Contribution (in %) of different background components in the sideband, signal region 
and for the wrong-charge combination for the (717, ir s ) pair. 



component 


sideband 


signal 


wrong charge 


qq, q = u, d, s 


40.9 


17.9 


40.8 


cc 


17.4 


30.6 


19.1 


B + B~ generic 


24.3 


28.2 


21.3 


B°B° generic 


17.4 


23.3 


18.8 



For the lifetime measurement, the combinatorial background parameters have been fitted to 
the data for wrong-charge combinations with M m i ss > 1.84GeV/c 2 . As a cross-check, a sideband 
region (1.8 < M m i ss < 1.84GeV/c 2 ) is also defined. The signal region is taken to be the interval 
M m i ss > 1-854 GeV/c 2 . The fractions of the different components in the sideband, signal and wrong 
charge region are shown in Table [l]. 

3.2 Lifetime measurement 

3.2.1 B° — ► D*~ 7T+ vertex determination 

The B° — > D*~-k + vertex has been computed fitting the fast and slow pion tracks with the beam 
spot constraint. To take into account the flight of the B meson in the xy plane, the effective size 
of the luminous region in the vertical direction has been assumed to be 60 /im. 

The resolution on the B° — ► D*~ir + vertex is dominated by the beam spot and by the fast pion 
track. Therefore, we have required the fast pion to have at least 3 SVT hits. The longitudinal 
or z precision on this vertex is strongly correlated to the direction of the fast pion track and is 
best when this direction lies in the transverse plane, in which case the resolution is 50 jttm. The 
error on the z coordinate of this vertex is required to be less than 150 //m. The z resolution of 
the B° -> D*-ir + vertex can be fitted with the sum of two Gaussian distributions, where the core 
Gaussian has a width of 48 fim and contains 80% of the events. 

3.2.2 Rejection of D° tracks and B tag vertex determination 

Before fitting the vertex of the second B in the event, hereafter named Bt ag , it is necessary to reject 
tracks coming from the non-reconstructed D°. Including these tracks in the fit would result in a 
bias, pulling the position of the Btag vertex systematically closer to the B° — > D*~ir + vertex and 
thereby reducing the measured lifetime. Thanks to the high momentum of the D , most of its decay 
products will be boosted in a cone around the parent direction. To reject these, tracks are excluded 
that lie within the region in the 6tk ~P plane that is populated by D° decay particles, where Otk 
is the the angle between the direction of a given track and the reconstructed D° direction and p 
is the momentum of the track in the center of mass frame. This requirement is applied up to a 
maximum angle of 2 rad. It has an efficiency of 45% for the tracks from the Btag and 1.8% for the 
D° decay products. As a result, 98% of the accepted tracks come from the B tag - 

The Btag vertex has been fitted from the tracks selected as explained above, discarding the 
tracks with excessive contribution to x 2 . To achieve a high quality vertex it is required that the 
fit be performed using at least two tracks. Even after this treatment, some tracks from the decay 
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of charmed hadrons are still included in the B tag vertex fit. As these come from the charm decay 
vertex, which is generally at larger z values due to the charm hadron lifetime and the boost in the 
laboratory frame, this produces a bias in the mean value of the z coordinate of the order of 20 /im. 
The z resolution for the Btag vertex can be fit with the sum of two Gaussian distributions, where 
the core Gaussian has a width of 92 ^m and it contains 70% of the events. 

3.2.3 Vertex quality selection 

Further selection cuts are applied to improve the vertex quality. The number of SVT hits associated 
to the fast pion is required to be at least 4, the z error on the B° — > D*~tt + vertex fit is required 
to be less than 150 fim, the normalized x 2 of the B° — > D*~ir + vertex fit is required to be less than 
20, and at least two tracks are required to be used by the Bf ag vertex fit. The efficiency for these 
cuts estimated from Monte Carlo simulation to be 58%. The number of signal events is 1696 ± 107 
and the fraction of combinatorial background in the signal region aback is 27 ± 6%. 

3.2.4 B° lifetime fit 

Using the fit result for the two B vertices, it is possible to extract a longitudinal vertex separation, 
Az 

— z D*ir ~ z tag- At is then computed using the approximation At — Az/P'yc, where the boost 
/?7 of the T(4S) at PEP-II is known to be 0.56. In order to fit the At distribution, three regions 
were defined in the M m i ss spectrum: a signal region for candidates with M m i ss > 1.854 GeV/c 2 , a 
sideband region with 1.800 < M miss < 1.840 GeV/c 2 and a wrong charge region for candidates with 
both pions have the same charge and M m i ss > 1.840 GeV/c 2 . 

The wrong-charge distribution is used to define the shape of the combinatorial background 
under the signal. The At distribution of the wrong charge region is fit with the function: 

/^(At) = (1 - a B )(a u d s gauss(At,a uds ) + (1 - a uds ) gauss(At,a c )) 

+ob((1 - a bw )E G (At - b n ; T back ,a bn ) + a bw E G (At - b w ;T bac k,^bw)) (1) 

where Eoit; r, a) is the convolution of an exponential with lifetime r with a gaussian of width a. 

The parameters a bw , b n , a bn , b w , a bw for the background "resolution function" are fixed to the 
resolution function determined from Monte Carlo simulation for signal events. The background 
lifetime is fixed to the value, 1.55 ps, measured from the same distribution for BB events in the 
Monte Carlo sample. The widths u u ds an d er c of the two Gaussian distributions used to parametrise 
the continuum light quark {uds) and charm components are fixed to the values, 0.70 and 0.96 ps 
respectively, extracted from the Monte Carlo. «e and a u d s are the only parameters left free in the 
fit of the background At distribution. 

To take into account the fact that the fractions of uds and charm in the different samples are 
different (cf. table |l|), a u d s is fixed to the value 0.37, obtained from the Monte Carlo simulation in 
the signal region, for the final lifetime fit. 

The lifetime fit is performed using an event-by-event resolution estimate obtained by scaling the 
error given by the vertex fit for the two vertices. The Monte Carlo Az pull distribution has been 
used for the rescaling. This can be described by two Gaussian distributions with widths p n = 1.04 
and p w = 2.51, centered at s n = —0.22 and s w = —0.95 respectively. The fraction a w id e of the 
events in the wide Gaussian is 11%. 

Therefore, the probability density function for the unbinned maximum likelihood fit for event i 

is: 

f(Ati, e^) = (1 - a bac k)fs(Ati, Oi) + a back f b (Ati) (2) 
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where 



f 3 (Ati,ai) = (1 - a wi de)E G (Ati - s n ai;T B o,Sp n ai) + a wi d e EG(Ati - s w af,T B o,Sp w (Ti) (3) 

Here At, and a% are At and its error coming from the fit of the two B vertices and S is a global 
scale factor for the errors. 

The result of the unbinned maximum likelihood fit with one free parameter, the B° lifetime, is 
t b o = 1.51 ± 0.05 ps as shown in Fig |I[ 




Figure 3: Missing mass distribution in par- 
tially reconstructed B° — > D*~ir + events af- 
ter all cuts. The data are represented by 
points. The contribution of the qq final 
states from Monte Carlo is: uds (black), c 
(dense hatched), charged B (light hatched), 
neutral B (white). 



Figure 4: Distribution of At in ps. The data 
are represented by the points with error bars. 
The continuous line shows the result of the fit 
and the shaded area shows the contribution 
of the combinatorial background (27% of the 
sample) . 



3.2.5 Bias 

Using the same procedure, a lifetime fit has been performed on a signal Monte Carlo sample, giving 
t b o = 1.48 ± 0.025 ps. Removing the tracks from the D°, the bias is removed and the result, 
t b o = 1.52 ± 0.02 ps, is compatible with the generated lifetime t b o = 1.56 ps. The estimated bias is 
then 0.04 ps and this has been added to the fit result from data to obtain the quoted final lifetime 
value. 

Table Q shows the fit results obtained with different values for the cone cut, after applying the 
appropriate Monte Carlo correction for bias. The corrected lifetime is stable with respect to the 
variation of the cone cut. 

3.2.6 Systematics 

A preliminary evaluation of the systematic errors on the lifetime result has been performed. 
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Table 2: Results of the fit obtained with different values for the cone cut for a smaller data sample 
corresponding to 5.3 fb _1 after applying the Monte Carlo bias correction. All the times are in ps. 



Cut (rad) 


Events 


7~raw 


At mc < 


7~corr 


1 


2595 


1.42 ± 0.04 


0.12 ±0.02 


1.54 


1.5 


1839 


1.49 ±0.05 


0.06 ±0.02 


1.55 


2 


1287 


1.52 ±0.06 


0.04 ±0.02 


1.56 



Table 3: Sources of the systematic error and their impact on r. 



Source 


Variation 


At (ps) 




27 ± 6% 


0.039 


OiB 


49 ± 3% 


0.012 


Quds 


37 ± 20% 


0.006 


Tback 


1.55 ±0.05 ps 


0.010 


s 


1.00 ±0.13 


0.051 


bias 


0.04 ± 0.02 


0.020 


Total 




0.07 



The combinatorial background fraction a\, ac k has been varied by 6% leading to a variation 
of 0.039 ps in the fitted lifetime. The parametrisation of the combinatorial background has been 
varied yielding variations of 0.012 {olb BB fraction), 0.006 ( a u d s uds fraction) and 0.010 ps (r& ac fc 
background lifetime). 

This analysis relies on the resolution function fit from the Monte Carlo simulation and therefore 
is sensitive to possible discrepancies between simulation and data. The global scale factor S has 
been left free in the fit giving S = 1.13 ± 0.12. The corresponding variation of the lifetime of 
0.051 ps has been added to the systematic error. 

The total systematic uncertainty of 0.07 ps, is dominated by the uncertainty in the continuum 
background fraction and the uncertainty in the vertex resolution function (table |3|) . As these two 
contributions are related to variations estimated on the data themselves, further increases in the 
size of the data sample should reduce these systematic errors. 

As a further cross-check of the resolution in Az, a data-Monte Carlo comparison has been 
performed for the sideband and for the off resonance data. The agreement is good in both cases. 

4 B° to D* £ i>i analysis and results 

4.1 Event and track selection and sample composition 
4.1.1 BB event selection 

Hadronic events were selected by requiring that at least four charged tracks be present. To reduce 
the contamination from light quark and cc events, the second Fox- Wolfram moment R2 is required 
to be smaller than 0.5. 
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4.1.2 Track selection 

Only charged particles tracks were employed for this analysis. High momentum tracks were recon- 
structed by matching hits in the Silicon Vertex Tracker (SVT) with track elements in the Drift 
Chamber (DCH). Low momentum tracks do not hit enough wires in the DCH due to the bending 
induced by the magnetic field. They are however reconstructed by the SVT stand-alone pattern 
recognition, with more than 50% efficiency for momentum p > 70 MeV/ c. 

4.1.3 Particle identification 

Charged tracks are identified as electrons using information from the electromagnetic calorimeter 
(ratio of the energy released to the associated track momentum, transverse profile of the electromag- 
netic shower), from the energy loss in the DCH and from the Cerenkov angle as measured with the 
Detector of Internally Reflected Cherenkov light (DIRC). The efficiency for electron identification 
in the acceptance of the electromagnetic calorimeter is about 90%, with a hadron mis-identification 
probability < 1%. 

Muons were selected by requiring penetration through most of the Instrumented Flux Return 
(IFR). The selection requirements result in ~ 75% efficiency for muon identification with ~ 2% 
hadron mis-identification probability. The Cherenkov light information in the DIRC was then 
employed to further reject mis- identified kaons, by requiring the likelihood for the kaon hypothesis 
to be less than five percent. 

Kaons are used for B° flavour tagging. Among all charged particles with momentum p > 
500MeV/c, kaon candidates are selected on the basis of the energy loss in the DCH and in the 
SVT, as well as observed Cherenkov angle and number of photons in the DIRC. True kaons were 
selected with ~ 80 % efficiency, while ~ 1% of the pions were mis-identified as kaons. 

4.1.4 B° — > D*~£ + vi sample selection 

Only the charged pion and the lepton are required to reconstruct B° — > D*~i + u^ decays. The 
lepton candidate is required to have momentum in the T(4S) rest frame (k* ) in the range 1.4 < 
k* < 2.3GeV/c, while the tt* had to satisfy the requirement k* < 190MeV/c. 

Events were then selected by exploiting the kinematics of the decay. The B° four momentum 
is computed by assuming that the B° is produced at rest in the T(4S) decay, thereby neglecting 
the small boost of the B meson in this frame, where its momentum is about 300MeV/c. Due to the 
limited phase space available in the decay D*~ — > D°tt~, the tt* is emitted inside a fairly restricted 
cone centered about the D*~ direction in the laboratory frame. The D*~ four-momentum can 
therefore be computed by approximating its polar and azimuth angles with those for the tt* , and 
parametrising its momentum as a linear function of the tt* momentum: 



The D*~ direction is determined with an accuracy of about 15°, the error on the D*~ momentum 
is about 400 MeV/ c. The neutrino escapes detection, but its invariant mass can be computed from 
the B°, D*~ and l + four-momenta from the relation: 



Neglecting resolution effects, this quantity must be zero, while background events are spread over 
a wide range of M 2 . 



p{D*~) = a + (3-p(?r*) 



(4) 



M u 2 = (V(B°) - V(D*~) - V{t + )) 2 



(5) 
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The B° decay point is determined by intersecting the it* and £ + tracks, constrained to the 
beam-spot position, which was determined on a run- by-run basis using Bhabha events. Only events 
for which the probability of the vertex fit was greater than 1% were further considered. This cut 
rejected about two thirds of the combinatorial background, while still retaining ~ 80% of the signal. 

4.1.5 Sample composition 

The combinatorial background is determined from the data using those events for which the £ and 
7r* have same electric charge ("wrong-charge"). The contribution of the events from any decay 
process of the kind B — > D*~£ + vX ("resonant events") in the mass band {M v 2 > — 2GeV 2 /c 4 ) 
was determined by subtracting the wrong-charge from the right-charge sample, normalized in the 
sideband region -10 < Mj 2 < -5GeV 2 /c 4 . A small correction (~ 5%) was applied to the wrong- 
charge on-peak events to account for events of the kind B°^£, B°—fD*, which populates the right- 
and wrong-charge samples differently. The number of events obtained in the T(4S) resonance data 
sample is reported in the first line of Table ||. The second line shows the number of events obtained 
in the continuum data sample, multiplied by the ratio of the resonant and continuum luminosities: 
non BB events account for about 15% of the combinatorial background in the mass band and 5% 
of the sample obtained after subtraction of the combinatorial. 

Table 4: Data statistics. The number of continuum events is normalized by the ratio of integrated 
luminosities. 



Data Sample 


Total 


Resonant 


On-peak 
Off-Peak 


190080 
23900 


89360 
4160 



In the Monte Carlo simulation the difference between the right charge and the rescaled wrong 
charge events in the mass band is — 2.1 ± 1.6% No correction is applied to the data for this effect, 
and a relative systematic error of ±2% is assumed. The number of resonant events in the sample 
is therefore: 

Mres = 89360 ± 500 stat. ± 1800 syst. (6) 

Figure || shows the M 2 distribution for the electron data sample after subtraction of the off-peak 
data, before (left) and after (right) subtracting the combinatorial background. 

Several physical processes contribute to the D*£v final state. Contamination from B + is 
mostly due to the decay B + ^D**°£u where the orbital charm excited state, _D**° 5 decays into 
D*~tc + . The fraction of these events is 9±5%, where the error includes the systematic uncertainty 
due to the efficiency of the B — > D**£ + ui selection. The composition of the sample is reported in 
Table |5| All the processes originating from a B° decay are considered as signal for this analysis. 

4.2 tbo measurement and mixing 
4.2.1 Proper time determination 

The position of the B° — > D*~£ + U£ ("decay") vertex is reconstructed as described above from the 
£, tt* tracks and the beam spot. The decay point of the other B ("tag" vertex) is determined using 
the remaining tracks in the event, constrained to the decay vertex in the plane orthogonal to the 
beam axis (x,y). This constraint exploits the fact that in the x-y view the separation between 



15 




Figure 5: distribution for the right charge data, after subtraction of the off-peak events for 
electron candidates. Dots with error bar: data; histogram: Monte Carlo simulation. The right 
hand side plot is obtained after subtraction of the combinatorial background. 



Table 5: Sample Composition in the signal region (in %). 



B" -> D*-£ + u e 


36.6 


B°^D** 


2.7 


B°^D*-X c /t 


0.4 


B^D*-fake I 


0.8 


B+^D** 


4.5 


off-peak resonant 


2.4 


D*/l from different Bs 


4.0 


Other BB Combinatorial 


40.7 


continuum 


7.9 



the B° and the B° is small compared to the experimental resolution; it is applied to remove the 
tracks from charm or hyperon decay from the vertex fit, which would bias the Az reconstruction. 
To reduce the additional bias, peculiar to inclusive analyses, due to the un-resolved particles from 
the D° produced at the decay vertex, all the tracks lying within a one radian cone around the ir* 
are rejected from the tag vertex fit. In this way, about 60% of the D° tracks are removed, while 
still retaining ~ 70% of the B tag tracks. 

The B° lifetime is determined by measuring event by event the distance along z between the 
decay and the tag vertex, as already described in section 3.2.4. 

In the Monte Carlo simulation the Az resolution for the signal events is described by two 
Gaussian distributions: the narrow (wide) component containing about 65% (35%) of the events, 
with a width of 120 (370) /im and a small bias of about 20 (65) /mi. The pull distribution is 
parametrised by the sum of two Gaussian distributions, the narrow (wide) component contained 
82% (18%) of the events, with a width of 1.07 (2.4) and a bias of 0.19 (0.7). The same resolution 
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was assumed for B + events. 



4.2.2 Determination of t b o 

The B° lifetime is determined by means of an unbinned maximum likelihood fit, accounting for the 
event-by-event error determined by the vertex reconstruction algorithm. The fit function was the 
sum of the probability density function (pdf) for B°, B + and combinatorial background multiplied 
by the relative fractions for each type in the sample. 

The Az probability density function for B° events is described by the convolution of the proper 
time exponential with the pull resolution function described above. The same function is assumed 
for B + events, with a lifetime fixed to the present world average t b + = 1.65 ± 0.04 ps. 

The combinatorial background is described by the sum of a zero-lifetime and a non-zero lifetime 
component, both convoluted with a double Gaussian resolution function. All parameters (lifetime 
and fraction of non-zero lifetime events, widths, fractions and biases of the resolution function) 
were determined independently in the data and in the simulation by fitting wrong-charge events 
with the same program used to determine t b o . Good consistency is found in the data and in the 
simulation, except for the rms width of the resolution function, which is 1.06 times larger in data. 
This discrepancy is attributed to resolution effects. 

The method was first tested on the dedicated sample of simulated B° — > D*~£ + v^ events after 
applying the selection requirements described above. Due to the influence of the residual tracks from 
the un-reconstructed D°, the measured lifetime in Monte Carlo is smaller than the generated one. 
The resulting reconstruction bias At = 0.217 ± 0.010 ps for the lifetime is defined as the difference 
between the generated value of 1.560 ps and the fitted value T B o raw,MC = 1.343 ± 0.010 ps. 

The final lifetime result for the data t b o was obtained by adding the bias At determined 
previously in the Monte Carlo to the raw fit result T B o raw,Dt = 1.405 ± 0.016 ps. The lifetime 
obtained in this way is: 

t b o = 1.622 ± 0.016 ps (7) 
where the error is statistical only. Figure |6| shows a comparison between the data and the fit result. 



4.3 Systematic studies 

The study of the systematic effects is still preliminary. The following sources of error have been 
considered: 

• Sample Composition. The amount of B + is varied by ±50%; the fraction of combinatorial 



background is varied by ±2 % according to the study described in Section 4.1; the fractions of 
events from uncorrelated I D*~ production falling in the right- and wrong-charge samples are 
also varied within their allowed range, as computed from the measurements of the branching 
ratios B(B—>£X), B(B—>D*~ X), and of the B°B° mixing parameter Xd- 

Background Lifetime. The B + lifetime is varied within its error; the parameters of the 
pdf representing the combinatorial background is varied within their allowed range, properly 
accounting for correlations: this gives an error of ±0.02 ps. Alternatively, the background 
parameters determined from a fit to the data side band are used, and the difference in the 
result (±0.07 ps) is added in quadrature to the systematic error. 

Resolution Effects. The fit was repeated by scaling the resolution function in the data 
by the factor 1.06 as indicated by the fit to the background sample: the measured lifetime 
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Figure 6: Lifetime fit result on a linear (left) and a logarithmic scale (right). The inner dark area 
is the B + contribution, the lighter area is the combinatorial background, the line is the result of 
the fit. Dots with error bars represent the data. 



decreased by 0.06 ps. This difference is considered as the systematic error due to the detector 
resolution. If the fraction of events in the narrow Gaussian describing the resolution function 
were varied by ±5% the fit result would vary by ±0.04 ps. This is currently taken as a rough 
estimate of error for the pdf in data. The stability of the result versus some relevant variables 
(^,ps> ^psjJtypsj M v 2 ) was also studied, but no significant effect was found. 

• Analysis Bias. The sizeable correction for bias is determined from the Monte Carlo sim- 
ulation. In the data, the fit was repeated for different values of the tt* cone cut. The 
systematic error was computed as half the difference between the largest and the smallest 
values obtained. 

Table ^| shows the components of the estimated systematic error. 

Table 6: Contributions to the estimated systematic error. 



Source 


Variation 


a(r B o) ps 


B + fraction 
B^D*-,B^£ 
Comb. frac. 
Background pdf 
t(B+) 
Error Scale 
Data pdf 
At 


±50% 

±2% 

±0.04ps 
1.06 

±5% 


±0.005 
±0.002 
±0.014 

+0.070 
-0.020 

±0.005 
-0.060 
±0.040 
±0.033 


Total 




+0.091 
-0.085 
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5 Summary 



A large sample of B° — > D*~ir + and B° — > D*~l + vp decays has been obtained from the data 
sample collected by the BABAR detector at the PEP-II storage ring by means of partial D*~ 
reconstruction. Only the low energy ir~ from the decay D*~ — > ir^D was tagged and coupled to 
a high momentum opposite charge particle, either a pion or a lepton. 

In the -B — > D* - 7r + case, the missing mass recoiling against the system of the two pions was 
determined using the beam energy constraint. A clear excess of 10922 ± 235 events in the D° mass 
range was observed, and interpreted as the signal of B° — > D*~it + production. 

Events from the B° — > D*~£ + i>£ three body decay were selected by reconstructing the mass 
of the un-observed neutrino: for this purpose, the B° four momentum was computed neglecting 
its small boost from the T(4S) decay, and the D*~ four momentum was computed by properly 
parametrising the one of the soft pion. A signal of (89.4 ± 0.5stat. ± 1.8syst.) x 10 3 events was 
observed. 

The lifetime of the B° meson was determined from the distance, projected along the beam 
direction, between the tagged and the recoil B° meson in the event. The preliminary results are: 

t b o = 1.55 ± 0.05 ± 0.07 ps (D*ir), 
t b o = 1.62 ± 0.02 ± 0.09 ps {D*lv), 

were obtained by the two analyses. The two values are in good agreement, while the mutual 
correlation is small. 

Preliminary studies indicate sufficient tagging efficiency and flavor separation for a time depen- 
dent measurement of the B° mixing process. Clear evidence of a time dependent asymmetry is 
observed in the data, consistent with the expectation from the simulation. This will allow a precise 
determination of the B°B° mass difference Am^ in the future. 
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